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Abstract 

Despite the intense theoretical and experimental effort, an understanding of the superconducting 
pairing mechanism of the high-temperature superconductors leading to an unprecedented high 
transition temperature T c is still lacking. An additional puzzle is the unknown connection between 
the superconducting gap and the so-called pseudogap which is a central property of the most 
unusual normal state. Angle-resolved photoemission spectroscopy (ARPES) measurements have 
revealed a gap-like behavior on parts of the Fermi surface, leaving a non-gapped segment known as 
Fermi arc around the diagonal of the Brillouin zone. Two main interpretations of the origin of the 
pseudogap have been proposed: either the pseudogap is a precursor to superconductivity, or it arises 
from another order competing with superconductivity. Starting from the t-J model, in this paper 
we present a microscopic approach to investigate physical properties of the pseudogap phase in the 
framework of a novel renormalization scheme called PRM. This approach is based on a stepwise 
elimination of high-energy transitions using unitary transformations. We arrive at a renormalized 
'free' Hamiltonian for correlated electrons. The ARPES spectral function along the Fermi surface 
turns out to be in good agreement with experiment: We find well-defined excitation peaks around 
u = near the nodal direction, which become strongly suppressed around the antinodal point. 
The origin of the pseudogap can be traced back to a suppression of spectral weight from incoherent 
excitations in a small w-range around the Fermi energy. Therefore, both mentioned interpretations 
of the origin of the pseudogap can not be held. Instead, the pseudogap is an inherent property 
of the unusual normal state caused by incoherent excitations. In a subsequent paper, also the 
supercunducting phase at moderate hole doping will be discussed within the PRM approach 18 . 

PACS numbers: 71.10.Fd, 71.30.+h 
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I. INTRODUCTION 



Since the discovery of superconductivity in the cuprates 1 , enormous theoretical and ex- 
perimental effort has been made to investigate the superconducting pairing mechanism which 
leads to an unprecedented high transition temperature 2^2— An additional puzzle is the 
unknown connection between the superconducting gap of the superconducting phase and 
the so-called pseudogap which is a central property of the most unusual normal state of 
the cuprates. In particular, the pseudogap has been subject to intense debates. Studies us- 
ing angle resolved photoemssion spectroscopy (ARPES) have revealed several key features 
of the pseudogap in the cuprates by elucidating the detailed momentum and temperature 
dependence^—. It was found that the pseudogap opens on a part of the Fermi surface (FS) 
around the anti- nodal point, leaving a nongapped FS segment known as a Fermi arc around 
the nodal direction. The pseudogap also smoothly evolves with decreasing temperature into 
the SC gap and was, therefore, interpreted in favor of a "precursor pairing" scenario^,—,—. 
On the other hand, there are several experimental and theoretical reports which suggest a 
different origin for the pseudogap, such as caused by another order which competes with 
superconductivity 8 . Superconductivity is usually understood as an instability from a non- 
superconducting state. Therefore, often in theoretical investigations, the starting point was 
either the Fermi-liquid or the anti-ferromagnetic phase at large or low doping. In this paper, 
we take a different approach and only consider hole fillings, in which either a superconducting 
or a pseudogap phase is present. 

A generally accepted model for the cuprates is the t-J model which describes the electronic 
degrees of freedom in the copper-oxide planes for low energies. Alternatively, one could also 
start from a one-band Hubbard Hamiltonian as a minimal model. However, for low energy 
excitations, the latter model reduces to the t-J model so that both models are equivalent. 
As our theoretical approach, we use a recently developed projector-based renormalization 
method which is called PRM^. The approach is based on a stepwise elimination of high- 
energy transitions using unitary transformations. We thus arrive at a renormalized 'free' 
Hamiltonian for correlated electrons which can describe the pseudogap phase. The obtained 
ARPES spectral function along the Fermi surface is in good agreement with experiment: We 
find well-defined excitation peaks around uj = near the nodal direction which are strongly 
suppressed around the antinodal point. The origin of the pseudogap can be traced back 
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to a suppression of spectral weight of the incoherent excitations in a small cu-range around 
the Fermi energy. Therefore, the usual interpretations of the pseudogap origin can not be 
held. Instead, the pseudogap is an inherent property of the unusual normal state caused by 
incoherent excitations. 

First, after a short introduction of the model in Sec. II, it seems to be helpful, to start 
from a short outline of the basic ideas of our theoretical approach (PRM) in Sec. III. A review 
of this approach has been given elsewhere^. Then, in Sec. IV, the PRM will be applied to 
the t-J model in order to investigate the pseudogap phase at moderate hole doping. The 
final results will be discussed in Sec. V. In a subsequent paper, the supercunducting phase 
will also be discussed. 

II. MODEL 

A generally accepted model for the cuprates is the t-J model. In particular, in the 
antiferromagnetic phase at small doping, it has turned out that it can be used to describe the 
electronic degrees of freedom at low energies. We adopt the same model also for somewhat 
larger hole concentrations, outside the antiferromagnetic phase, where the superconducting 
and the pseudogap phases appear 

7-t = - Uj c\ a c ja -fj.y~] c\ a c ia + ^2 jys&j = : tu + Wj. (i) 

ij,a i<7 ij 

The model consists of a hopping term 7i t and an antiferromagnetic exchange Jtj. Here, tij 
stands for the hopping matrix elements between nearest (t) and next-nearest (t 1 ) neighbors. 
Jij is the exchange coupling and \x is the chemical potential. The quantities 

c\ a = c| CT (l - Ui-a), ha = C ia (l - Tli-a) (2) 

are Hubbard creation and annihilation operators. They enter the model, since doubly occu- 
pancies of local sites are strictly forbidden due to the presence of strong electronic correla- 
tions. Note that the Hubbard operators restrict the unitary space to states with only either 
empty or singly occupied local sites. They obey nontrivial anti-commutation relations 

[c\ a ,c ja/ ] + = S ij (8 aa/ 'D a (i) + Sa-ciS?), (3) 

where the operator 

V c {i) = 1-v (4) 
3 



can be interpreted as a projector which projects on the local subspace at site i consisting 
of either an empty or a singly occupied state with spin a. Finally, n ia = <y ia Ci a is the local 
occupation number operator for spin a, and Sf is the a — ±1 component of the local spin 
operator 

Si = \^^pc\ a c m (5) 

a/3 

where a a p = a^p e f ls ^ ne vector formed by the Pauli spin matrices. In Fourier notation, 
the t-J model (pQ) reads 

H = 5^^ k " V) cLc k(7 + ^ ^qS q S_ q = H t + Hj, (6) 
k.cr q 

F , - _ f.. P M*4-nj) i - ST r.. P »q(Ri-B.j) 



Note that for convenience, we shall somewhat change the notation. From now on, all energies 
will be measured from the chemical potential, i.e., — /i will be denoted by e^. 



III. PROJECTOR-BASED RENORMALIZATION METHOD (PRM) 

Let us start with a short introduction to the projector-based renormalization method 
(PRM)^ 1 ^ which we shall use as our theoretical tool. The general idea is as follows: The 
method starts from a decomposition of a given many-particle Hamiltonian 

n = n + n 1 (7) 

into an unperturbed part Tio and a perturbation T-i\. In TCi, no parts should be contained 
which commute with TIq. Therefore, Tli accounts for all transitions with non-zero energies 
between the eigenstates of Tio- The aim of the PRM is to construct an effective Hamiltonian 
which has the same eigenspectrum as 7i, and which can be solved. The first step is to 
construct a new renormalized Hamiltonian 1~C\ which depends on a given cutoff A, 

H\ = Ho : \ + Hi,x, (8) 

with renormalized parts Tlo,\ an d 1~ti,\- Thereby, 7i\ should have the following properties: 
(i) The eigenvalue problem of 7Y ,a can be solved 

H ,x\n x ) = E x \n x ), 
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where E x and |n A ) are the renormalized eigenenergies and eigenvectors, (ii) From TCi^x, all 
transition operators are eliminated which have transition energies (with respect to 7^o,a) 
larger than the cutoff energy A. As shown in Refs.— the renormalization step from 7i to 
1~L\ can be done by use of a unitary transformation. Therefore, the eigenspectrum of 7i\ is 
the same as that of 7i. 

The realization of the renormalization starts from the construction of 7i\. Here, the 
knowledge of the eigenvalue problem of 7io,A is crucial. It can be used to define generalized 
projection operators, Pa and Qa, 

P A -4 = ^|n A )(m A |(n A |^|m A )e(A-| J E n A -^|), 

771,71 

Q A *4 = (1 - P X )A, (9) 

which act on usual operators A of the Hilbert space. Note that in Eq. (Q the vectors \n x ) 
and |m A ) are necessarily neither low- nor high energy eigenstates of 7^o,a- Pa projects on the 
part of A which consists of transition operators |n A )(m A | with excitation energies \E X — E^\ 
smaller than A, whereas Qa projects on the high-energy transition operators of A. 

In terms of Pa and Qa, the property of 7i\, not to allow transitions between eigenstates 
of 7^o, a with energy differences larger than A, reads 

q a ^a = o or n x = p x n x . (io) 

The effective Hamiltonian 7i\ is obtained from the original Hamiltonian 7i by use of a 
unitary transformation, 

H x = e x *He- x \ (11) 

where Xa is the generator of the unitary transformation, and the condition (flUj) has to be 
fulfilled. The renormalization procedure starts from the cutoff energy A = A of the original 
model 7i and proceeds in steps of width AA to lower values of A. Every renormalization 
step is performed by means of a new unitary transformation, 

^(a-aa) = e x ^H x e' x ^\ (12) 

Here, the generator X Aj aa of the transformation from cutoff A to the reduced cutoff (A — AA) 
has to be chosen appropriately (see below). In this way, difference equations are de- 
rived which connect the parameters of H\ with those of 7Y(a-aa)- They will be called 
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renormalization equations. The limit A — > provides the desired effective Hamiltonian 
7i = 'Hx^q = 7io,A^o- The elimination of all transitions in the original perturbation Tli 
leads to renormalized parameters in 7i ,A^o- Note that H is diagonal or at least quasi- 
diagonal and allows to evaluate all relevant physical quantities. The final expression for Ti. 
depends on the parameter values of the original Hamiltonian TC. Note that 7i and 7i have, 
in principle, the same eigenspectrum because both Hamiltonians are connected by a unitary 
transformation. 

What is left, is to find an appropriate expression for the generator X AiAA of the unitary 
transformation which connects H.\ with 7Y(a-aa)- According to Eq. ( flOl) . X Aj aa is fixed by 
the condition Qa-aa^a-aa = 0. As is shown in Refs.— one can find a perturbation 
expansion for X AjAA hi terms of Tii. The lowest non- vanishing order reads 

= A [Q(a-aa)Wi,a] + " • • • (13) 

L 0,A 

Here, L 0j a is the Liouville operator, defined by the commutator \- ^A = [7^o,a,^4], for any 
operator quantity A. Note that Eq. ffl3|) can further be evaluated, in case the decomposition 
of Q(a-aa)7^i,a into eigenmodes of L 0j a is known. Formally written, we decompose 

Q ( a-aa)Wi,a = 5>5U A) where L , A F^ AA = ^ AA F^ AA , (14) 

V 

so that A^ AA is given by 

v A,AA 



IV. APPLICATION TO THE t-J MODEL 
A. Renormalization ansatz 

Our aim is to apply the PRM to the t-J model which is a generally accepted model 
for the low-energy properties of the cuprate superconductors. We consider a regime with 
moderate hole-dopings. The hole concentrations should be large enough for the system to 
be outside the antiferromagnetic phase but small enough to be in the metallic phase. Our 
first aim is to find the decomposition of the Hamiltonian into an 'unperturbed' part Tio and 
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into a 'perturbation' Tii. We assume that the hopping element t between nearest neighbors 
is large compared to the exchange coupling J. Therefore, Tit is the dominant part of the 
Hamiltonian in the metallic phase and should be included in 7i . However, also Tij has a 
part, which commutes with the hopping term, and which will be called Ti.^. Note that this 
part of Tij will not lead to transitions between the eigenstates of 7i t - Therefore, 7i t and 
Tij together form the unperturbed Hamiltonian 7Yo- The remaining part of TCj does not 
commute with Ti, t and forms the perturbation 7i\. Thus, we can write 

U -IV i 1-1 1-1 

In the framework of the PRM, the perturbation Ji\ will be integrated out by use of a 
unitary transformation. In lowest order perturbation theory, the generator of the unitary 
transformation X\&\ is given by Eq. ffTol) and relies on the decomposition of rij into the 
eigenmodes of L . However, it will be impossible to find the exact decomposition of TCj, due 
to the presence of Hubbard operators in 7i t - Therefore, we have to apply approximations. 
For this purpose, we start by decomposing the electronic spin operator 

* a/3 i 

into eigenmodes of L t instead of into eigenmodes of l_o. Here, U is the Liouville operator 
corresponding to the hopping part 7i t of H-q. The exchange Tij is given by a sum over 
products of spin operators S q ■ S_ q . Therefore, the decomposition of S q into eigenmodes of 
L t can be used to find an equivalent decomposition of Hj. 

The easiest way to decompose S q is to derive an equation of motion for the time-dependent 
operator S q (t), where the time dependence is governed by Tit, 

S q (t) = e Htt S q e~ iHtt = e iLtt S q . (17) 

Due to Eq. (E]), the first time derivative reads 

j t % = -^EtS^ 6 ^-^ (18) 



a/3 ij^l 

1 E%E*« eiqRi ( 1 - eiq(Ri " Ri) )4^. 



op Ijtl 

It can be interpreted as the hopping of a hole from some site / to a neighboring site % and 
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vice versa. The second derivative is characterized by a twofold hopping, 

4 (e^ - e***) (S^ (i) - S^o(0) (19) 



d 2 



^?EE E fa%(e'"W< lR ' 



2v ^ a/3 i? y i( ^) 

x {a aP (c) a V a {l) c ip + c)_ a Sfcipj + Ka {fipVail) c ja + c ] if3 S[ 

It has two different contributions. The first one describes the hopping of the hole from i 
back to site I from which it originally came and, equivalently, the hopping from I back to i. 
The second term in Eq. ( I19p stands for a twofold hopping away from the starting site. 
Let us discuss the first contribution to Eq. f[T9l in more detail. The operators 

V {i) = (l-n i ,i)0--n i ,i) (20) 

and Vq(1) can be interpreted as local projectors on the empty state at site i and site I, 
respectively. They assure that the original sites % and I were empty before the first hop. 
Their presence results from the fact that doubly occupancies of local sites are strictly for- 
bidden which is a consequence of the strong correlations in the t-J model. In a further 
approximation, let us replace Vo(i) and Vo(l) by their expectation values, 

V (i) => ((1 - ^, T )(1 - n u )) =: P , (21) 

which can be interpreted as the probability for a local site to be empty. Without the second 
term in Eq. ffl9l) . we are led to the following equation of motion for S q (i): 

|,S q = -^S, (22) 

Obviously, the differential equation (T221 describes an oscillatory motion of S q (£) with fre- 
quency u; q , where 



-I 



2Po(t q=0 " tl) = > 0, t\ = 4 e*<*~**\ (23) 



Note that the averaged projector P = 1 — n also agrees with the hole concentration 5 away 
from half-filling, i.e. Po = 5 = 1 — n, where n is the electron filling. 

Before carrying on with the physical implications of Eqs. (j22|) . (|23|) . let us discuss the 
influence of the hole (or electron) hopping in Eq. ffl9l) to second nearest neighbors and also 
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to more distant sites. As long as the dynamics of S q (t) is alone governed by the hopping 
Hamiltonian Tit, all these hopping processes are important and would have to be taken into 
account. For instance, for a state close to half-filling outside the antiferromagnetic regime, a 
hole and a neighboring electron can freely interchange their positions for a system governed 
alone by 7i t . The hole can easily move through the lattice. However, the situation is different 
from the case, for which the dynamics is governed by Tio = Tit + Ti^p . Then, we have to 
decompose the perturbation TCi into eigenstates of L , where L is the Liouville operator 
corresponding to Ho. Thus, the dynamics of S q is not governed alone by the hopping 
Hamiltonian TC t but also by the yet unknown commuting part H^ of Hj. However, in 
Appendix [A] it is shown that local antiferromagnetic spin fluctuations due to H^ restrict 
the hole motion to neighboring sites. The hopping to more distant sites is strongly suppressed 
by spin fluctuations. Therefore, the former equation of motion (122]) for S q (£) turns out to be 
a good approximation for the case that the dynamics is determined by the full unperturbed 
Hamiltonian Hq including the exchange part. 

The arguments in Appendix |A] are based on the evaluation of the dynamical spin sus- 
ceptibility x(q, a;) as follows. Using the Mori-Zwanzig projection formalism x(q, u) can be 
written as 

— u 2 

Here, uj^ ~ is approximately the frequency, given in Eq. (|23|) . and £ q (w) is the selfenergy. 
The exact expression of £ q (cu) in terms of the Mori scalar product reads 

S q M = . -. (QSJ -J— QS q ). (25) 

ql ; (S q |S q ) 1 ql cu-QL Q-*77 q; 1 ; 

Here, Q is a generalized projection operator which projects perpendicular to S q and S q 
(for details see Appendix D). Due to construction, the operator QS q in the 'bra' and 'kef 
of Eq. (125]) corresponds to the second line in Eq. (fl9j) . and describes a twofold hopping 
away from the original site. Therefore, the selfenergy £ q (u;) provides information about 
the hopping processes between next nearest neighbor sites and to more distant sites. In 
Appendix [A] the selfenergy S q (u;) is evaluated in a factorization approximation by including 
the spin fluctuations from Ti.^. The result is shown in Fig. CD, where the imaginary part 
of £ q (u;) for a small q- vector is plotted (solid line) in the presence of spin fluctuations due 
to Hj . As is seen, £ q (co>) is rather small and almost ^-independent over a wide frequency 




FIG. 1: Imaginary part of the selfenergy 9£ q (iij) from Eq. (|25p in the presence of spin fluctuations 
(J = 0.2t, solid line) and in the absence of spin fluctuations (J = 0, dashed line). The q- vector is 
fixed to q= (vr/20, tt/20). 

range. Thus, the only effect of S q (o;) is to give rise to a small damping and lineshift of the 
resonances of x(q, We have also repeated the same calculation for 9£ q (cc;) in the absence 
of 7ij\ i.e. when Ho is replaced by Tit (dashed line in Fig. CD). A strong tu-dependence is 
found for small q-values around u = 0. This shows that long reaching hopping processes 
are important in this case. From these findings, one can conclude that the hopping to 
more distant than nearest neighbors is of minor importance as long as the exchange part 
7ij ^ is not neglected in H, . A possible explanation would be that local antiferromagnetic 
correlations are still present at moderate hole doping outside the antiferromagnetic phase. 
They lead locally to strings of spin defects which are well known from the hole motion in 
the antiferromagnetic phase. 

Let us come back to the discussion of the oscillation behavior in Eq. ( |T9i) which can 
be understood as follows. When an electron hops to a neighboring site, it preferably hops 
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back to the original site, since this was definitely empty after the first hop. In contrast, the 
hopping to next nearest neighbor sites is energetically unfavorable due to local antiferromag- 
netic order. As will be shown in a forthcoming paper—, the proportionality of lj^ ~ 5 turns 
out to be the basic feature for the understanding of the superconducting pairing mechanism 
in the cuprates. The oscillation becomes less important for larger S which agrees with the 
weakening of the superconducting phase for larger hole doping. 
The solution of Eq. (1221) is easily found, 

S q (t) = S q COS UJqt + — S q SU1 Cj^t (26) 

= ^(S q - - S q ) + i(S q + - S q ) e-^\ 

where S q = S q (t = 0) and S q = |S q (t = 0) was used. From Eq. fl26l) . the decomposition of 
S q into eigenmodes of L can immediately be identified, 

Lo [ks q T -f S q )] = ±^ q [i(S q =F -f S q )], (27) 

Z UUn U^q 

which leads to the intended decomposition of the exchange TCj as follows: 

Hj = J n s q s_ q = J * (M<i) + M<i) + 4(q)) , (28) 

q q 

where 

Ml) = \ (s q S-„ + -^S q S_ q ) , (29) 

and 

LoA(q) = 0, LoA(q) = 2^ q A(q), L 4(q) = -2£ q 4(q). (30) 

Here, an additional approximation was used. In deriving Eqs. (130]) . the eigenmodes of the 
two spin operators S q • S__ q in the expression for Tij were taken separately from Eq. (|27|) . In 
this way, all local configurations were disregarded, where two spin operators in local space 
are located on neighboring sites. Thereby, a possible hopping between the two sites would be 
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obstructed. The inclusion of these processes would need additional considerations. However, 
they would not change our results substantially. 

With Eqs. (1251) . we have arrived at the intended decomposition of the t-J model. The 
Hamiltonian 

H = yi g kCA°- + y] JqSqS_ q (31) 
kcr q 

can be decomposed into an 'unperturbed' part Ho and into a 'perturbation' H±. It reads 
H = H t + H ,j =: 5^£k4<A<x + </ q A(q), 

kcr q 

Hi = (A(q)+4(q))- (32) 

q 

The aim of the projector-based renormalization method (PRM) is to eliminate all tran- 
sitions between the eigenstates of Ho which are induced by H\. Let us assume that all 
excitations with energies larger than a given cutoff A have already been eliminated. Then, 
the renormalized Hamiltonian H\ should have the form 

H\ = ^ £k,A 4<A(T + Jq,X PASqS_q , (33) 
kcr q 

however, with A-dependent prefactors e^x and J q A- Moreover, a projector was intro- 
duced which acts on operator variables. It guarantees that only transitions with excitation 
energies smaller than A remain from S q S_ q . 

The separation of H\ into an unperturped part 7i ,A and a perturbation Hi r \ reads in 
analogy to Eq. ([32]), H\ = H ,x + Hi : \, with 

Ho,\ = Ht,\ + Jq,A A),A(q) + E\, 
q 

n hX = J q.A e(A - |2< A |) (j^m + 4,A(q)) > ( 34 ) 

q 

where we have used the A-dependent extension of relation (l3"Uj) in order to exploit the prop- 
erties of Pa- Note that the 6-function 6(A — |2u) qi A|) in H\ t \ guarantees that only excitations 
with transition energies |2u> q A | smaller than A contribute to H\,\. In Eq. ( |3^|) . H t .\ is the 
renormalized hopping term from Eq. (1331) . H t< \ = J2ka £ k,A c^cW Also, the parameters 
<^q,A5 ^q,A and E\ in Eqs. ( l34l) now depend on A. Moreover, the new operators *4. Q; A(q) 
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(a = 0, ±1) depend on A, 



A,A(q) = j ( S q - ^- *J\ ( S_ q - ±- S_ q , A ) , (35) 

4 lA (q) = t f s q + -f- a J) ( s_ q + -J- s_ q , A ) , 

where u> qjA and S qjA are defined by 

^ = 2Po (^o,A - C)' C = E 4,a e^"^, 

*(#) 

S q ,A = ^[Wq,A, S q ] « -[W t) A, S q ]. (36) 

B. Generator of the unitary transformation 

To derive renormalization equations for the parameters of 7i A , we have to apply the 
unitary transformation (fl2l) to 7i\ in order to eliminate excitations within a new energy 
shell between A and A — AA. We use the lowest order expression ffTol) for the new generator 

x a ,aa = £ ^-e q (A, AA) (Xa(q) - 4,A(q)) • (37) 

q q ' A 

Here, q (A, AA) denotes a product of two 6-functions, 

6 q (A,AA) = e(A-|2< A |)e(|2a; qiA _ AA |-(A-AA)), 

which confines the elimination range to excitations with |2cj qjA _A A | larger than A — A A and 
|2c2) qiA | smaller than A. Roughly speaking, for the case of a weak A-dependence of |co> q;A |, 
the elimination is restricted to all transitions within an energy shell between A — AA and A. 
With (133]) . the generator X A aa can also be expressed by 

X XA x = 7^ q( A > AA ) (X S-q,A + S q , A S_ q ) . (38) 

q q>A 

In the following, we restrict ourselves to the lowest order renormalization processes. Then, 
J q A will not be renormalized by higher orders in J, and we can use J q A = «/ q from the 
beginning. 
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C. Renormalization equations 



The unitary transformation ([121) . applied to the renormalization step between A and 
A — AA, will be evaluated in perturbation theory in second order in J q , 



-1/(2) , 
^A-AA + 



(39) 



where 



n 

H 

n 



(0) 

A-AA 
(1) 

A-AA 



(2) 

A-AA 



S ) J ^k,A c^Cko- + E\ = H t ,x + Ex, 

kcr 

J q A,A(q) + [X XA x,H ttX \ + J q ( A " I 2 <a|) (A,A(q) + ^U(q)) 
q q 

2^A,AA, [Aa,AA,7^,a] ] + ^ -/q [^A,AA; " / ^0,a( { 1)] 

q 

+ E J q ( A - l 2 ^|) [ X ^.aa> A,A(q) + 4, A (q) ]• 



(40) 



Let us first evaluate Ti-^-Ax f rom second order processes. The commutators in Eq. (14"U|) 
are explicitly evaluated in Appendix A. Then, we can compare the obtained result with 
the formal expression for Tix-Ax which has the same operator structure as Tix, with A is 
replaced by A — AA. One obtains the following renormalization equation from the second 



order contributions in J q : 



£k,A-AA — £k,A 



^ ^2 ®q( A ' ( £ k+q,A + ^k-q,A ~ 2 £k,A) (S q ,A S_ qi A) 

q q>A 



X 



^£(i|) 2e <AAA)( £k ,- £k - q , f 

^2i^k',X - £k'+q,A - £k'-q,A)(c k ' -'C k / (7 ') 



(41) 



n 



(NL) 
k— q« ' 



where we have defined 



n- 



(NL) 
k,cr 



(4A*> - ^5Z(4 CT Ckv) 



(42) 



as non-local part of the one-particle occupation number per spin direction. An equivalent 
equation also exists for Ex- ax- Note that in Eq. fj4"Tj) an additional factorization approxima- 
tion was used in order to extract all terms which have the same operator structure as Tix- 
The quantity (S q AS_ q A) is a correlation function of the time derivatives of S q which can 
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easily be evaluated from Eq. (1B3I) . Note that an additional contribution to £k,A-AA, propor- 
tional to the correlation function (S q • S_ q ), has been neglected. The remaining expectation 
values in Eq. (T4"Ti) have to be calculated separately. In principle, they should be defined with 
the A-dependent Hamiltonian because the factorization approximation was employed for 
the renormalization step from Tix to 7^a-aa- However, 7i\ still contains interactions which 
prevent a straight evaluation of A-dependent expectation values. The best way to circumvent 
this difficulty is to calculate the expectation values with the full Hamiltonian Ti instead of 
with 7i\. In this case, the renormalization equations can be solved self-consistently, as will 
be discussed below. 

Note that the renormalization (J3TJ) of £k.A wa s evaluated from the second order part 
^a-aa °f the Hamiltonian (14"U1) . Thus, we are led to 

T~L\-a\ = 7~tt,\—A\ + 7~L\-a\ + -^a-aa, (43) 

where TCt,x-AX = J2k a £ k,A-AA cj^Cko-. What remains is to evaluate the renormalization part 
Hx-ax m fi rs t order in J q to 7Ya-aa- First, the second term on the right hand side of Eq. (140]) 
can be rewritten, since 

[X XA x,n t , x ] = -^J q q (A,AA) (X A (q) + 4 A (q)) . 
q 

Then, by combining the second and third term, we find 

^?-AA = ^J q A,A(q) (44) 
q 

+ J q ( A " l 2 ^|) e(A - AA - |2c£; q , A _AA|) (A,A(q) + 4^)) • 
q 

The excitation energies of Ai^q) and *4.j A (q) are restricted to |2u) q)A | < A by the first 
G-function in Eq. f j44l) . This condition is automatically fulfilled by the second 0-function, 
in the case that |2u> qi A-AA| only weakly depends on A and we can replace A by A — A A. By 
introducing the projector Pa-aa on all low-energy transition operators with energies smaller 
than A — AA, we find 

«?2aa = E J i P ^aa + A,A(q) + 4,x(q)) 

q 

= E^-^VS-q, (45) 
q 
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where we have used the representation (1281) for the scalar product S q • S_ q , 

S q • S_ q = A,A(q) + A,A(q) + A\ )X (q). (46) 
Finally, for the total Hamiltonian ^a-aa; we obtain according to (j4"3j) 

7~L\-A\ = ^ ^k,A-AA C^Ckcr + 2^ ^ q PA-AA S q • S_ q + E\—A\- (47) 
k,cr q 

Note that this expression completely agrees with the Hamiltonian at cutoff A, when A is 
replaced by A — AA. The required decomposition into Ho,a-aa and H^a-aa is found as 
follows. We use again the relation (f4"6l) . with A is replaced by A — A A, and rewrite Hx-ax as 

^i-AA = J 1 P ^ A * (•A).A-AA(q) + A,A-AA(q) + 4,A-Aa(<J)) • ( 48 ) 

q 

Using again Eq. fj45|) . we arrive at the renormalized Hamiltonian Hx-ax = 7^o,a-aa+7^i,a-aa 
in the following form, 

'Hq.x-ax = T~~£t,X—AX + 2_s ^q^o.A-AA(q) + E\- AX, 

q 

^1,A-AA = J( l @ ( A " AA " Ka-Aa|) (^l,A-AA(q) + 4,A-Aa(o)) ■ ( 49 ) 

q 

As expected, the renormalized Hamiltonians 7^ ,a-aa and Hi^x-ax have the same operator 
structure as at cutoff A. Therefore, we can formulate a renormalization scheme as follows: 
We start from the original t-J model, where the energy cutoff is denoted by A = A. Starting 
from a guess for the unknown expectation values, which enter the renormalization equation 
f|4"T]) . we proceed by eliminating all excitations in steps AA from A = A down to A = 0. 
Thereby, the parameters of the Hamiltonian change in steps according to the renormalization 
equation (14"TT) . In this way, we obtain the following model at A = 0: 

H\=o = H t ,x=o + 2_j ^qPA=oSq • S_ q + -£^=0 (50) 
q 

= £k,A=o cL c kCT + ^2 J q A),A=o(q) + E x =o- 

ktr q 

Note that in Eq. (loT)]) the perturbation Tii is completely integrated out. Only the part of 
the exchange, which commutes with the hopping term, remains. 

Unfortunately, due to the presence of the Ao-teim., the Hamiltonian 7ix=o can not be 
diagonalized. It does not yet allow us to recalculate the expectation values. Therefore, a 
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further approximation is necessary which consists of a factorization of the second term 

1 



q q V 



SqS_q + — S q .A=oS- q> A=0 ) • (51) 



According to Appendix B, 7Ya=o can finally be replaced by a modified Hamiltonian which 
will be denoted by 

n w = E^4^ + ET s q s - q + ^. (52) 

k(j q 

where the electron energy is modified according to 



- T7 ^2 A.",2 q ( £ k,A=0 - ^k-q,A=o) 2 ™ k -2<7> ( 53 ) 



iV ^ 4c^ A=0 



and is defined in Eq. (|42p . Note that the operator structure of agrees with that of 
the original t-J model of Eq. fl3Tj) . However, the parameters have changed. Most important, 
the strength of the exchange coupling in Eq. (|52|) is decreased by a factor 1/2. This property 
allows us to start the whole renormalization procedure again. We consider the modified t-J 
model of Eq. (|52|) as our new initial Hamiltonian, which has to be renormalized again. The 
initial values of Ji^ at cutoff A = A are £ k and J q /2. After the new renormalization 
cycle the exchange coupling of the new renormalized Hamiltonian Tv 2 ^ is again decreased 
by a factor 1/2, till after a sufficiently large number of renormalization cycles (n — > oo) the 
exchange operator completely disappears. Thus, we finally arrive at a 'free' model 

H = J2^c{ a c ka + E, (54) 

ko- 

where we have introduced as new notations = fi^ 00 ) ) 5 k = e£^ ca \ and E = f?( n->0 °). 
Note that the Hamiltonian 7i now allows us to recalculate the unknown expectation values. 
With the new values, the whole renormalization procedure can be started again till, after 
a sufficiently large number of such overall cycles, the expectation values have converged. 
The renormalization equations are solved self-consistently. However, note that the fully 
renormalized Hamiltonian (!54j) is actually not a 'free' model. Instead, it is still subject to 
strong electronic correlations which are built in by the presence of the Hubbard operators. 
Therefore, to evaluate the expectation values, further approximations have to be made. 



17 



D. Evaluation of expectation values 



The expectation values in Eqs. (I4ip and (|53|) are formed with the full Hamiltonian. To 
evaluate expectation values for operator variables A, we have to apply the unitary transfor- 
mation also on A, 

(A) = ^ e _ m ) = (A(X)) Hx = (A)„ , (55) 

where we have defined A(\) = e Xx Ae~ Xx and A = A(X — > 0). Thus, additional renormal- 
ization equations for *4(A) have to be derived. 

As an example, let us consider the angle-resolved photoemission (ARPES) spectral func- 
tion. It is defined by 

i r°° 

A(k, u) = — \ (cl(-t) c k(7 ) e^dt = <cL 8(1 + cj) c kCT > (56) 

Z7r J -oo 

and can be rewritten by use of the dissipation-fluctuation theorem as 

A(k,u) = ^^-csGfccj) , (57) 
where ^(^(k, uo) is the dissipative part of the anti-commutator Green function 

9G(k,a;) = — / (icU-t), c ka ] + )e^dt = ([cL,5(L + o;)c kCT ] + ). 

^ J -oo 

The time dependence and the expectation value are formed with the full Hamiltonian 7i, and 
L is the Liouville operator corresponding to 7i. According to Eq. (1531) . the anti-commutator 
Green function can be expressed by 

9fG(k,u;) = <[ci CCT (A),5(L A + ^)c kCT (A)] + > A , (58) 

where now the creation and annihilation operators are also subject to the unitary transfor- 
mation. To evaluate A(k,u), we have to derive renormalization equations for c k(J (A) and 
c ko .(A). According to Appendix [C] the following ansatz for c kcr (A) can be used: 

C kCT (A) = M kiA C k(T + — ^ V k; g ; A ^2 ^Z^ Q/3 ' ^WX^A ~ £k'+q,A) cL+qa C k '/3 C k+q7 . 

qk' 1' A a/?7 

(59) 

It can be justified from lowest order perturbation theory. Note that the A-dependence is 
transferred to the parameters u^x and f k , q ,A- Also the quantities u^x and e^x depend on 
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A. However, having in mind perturbation theory in J, this A-dependence will be neglected 
in the numerical evaluation of Sec. |V] below. According to Appendix [Cj the renormalization 
equations for w kA and fk, q ,A read 

ul 

x {n k / +q m k / + n k+q (D + n k ' - n k / +q )} 



4.A-AA = M k,A - ^2 (4^) q( A ' AA )( £ k',A -£k'+q,A) 2 j(^) + WM^aJ 



+ 4AT2 Yl 4^2 4^~( £k+( l'' X ~ £ k+q+q',A) (£k+q,A ~ ^k+q+q',*) 



qq 



q q' 



and 



x {v^, x Q q (\, AA) + ^ k , q , A e q ,(A, AA)} ^ 

X W+q'(™k+q+q' ~ n k+q — D) — m k+q n k+q+q / } (60) 



^k,q,A-AA = ^k,q,A +Mk,A@q(A, AA). (61) 



The quantities n k and m k in Eq. fl60|) are the k-dependent occupation numbers for electrons 
and holes per spin direction, which are formed with the full Hamiltonian 7i, 

n k = (4aCka), "Ik = (c kCT C kCT )- ( 62 ) 



In the following, we simplify the notation by suppressing the spin index a in (162]) . The 
renormalization equations (I60l) and (IBT!) for w kA and t>k,q,A, together with the ansatz (|59|) 
for c k)(J (A), enable us to evaluate n k and m k and also the ARPES spectral function. With 
some initial guess for n k and m k , we start from the parameter values of the original model 
at A = A, 

Mk,A = 1 , ^k,q,A = , (63) 

and eliminate all excitations in steps A A from A = A to A = 0. We end up with renormalized 
parameters which obey 



«k,A=0 7^ 1 > ^k,q,A=0 7^ °. 



Thus, after the renormalization, the annihilation operator c k (A = 0) =: cj^ at A = has the 
final form 

C k<r = M k,A=0Ck CT + 7^T7 2^ W k,q,A=0 "TTTJ / ^off ' 0V 7 J (£k',A=0 ~ £k'+q,A=0j C k/+qa C k '/3Ck +q7 - 



2N ^ ^ 

qk' q Q/37 
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As was discussed before, the Hamiltonian after the first renormalization Tv* 1 * can not directly 
be used to recalculate the expectation values n k and m^. In there is still a part of the 
exchange present, which is, however, reduced by a factor 1/2. Therefore, the renormalization 
has to be done again by starting from as the new initial Hamiltonian. Similarly. 



,(i) 

> u kcr 



can be considered as the new initial annihilation operator, i.e., c£2 = c k ^(A = A), with 

(i) (i) 

%,A=A - M k,A=0, %,q,A=A = W k,q,A=0- 

After n renormalization cycles, the exchange is scaled down by a factor (l/2) n . For the 
renormalization equation for and f k n qA , we obtain 

(41aa) 2 = (4:i) 2 -^E(i^) 2@ q ( A ' AA )(^-^) 2 (64) 



qk' V 1' 

X i l ) + 2^ V ^' X \ ^ k ' +C i mk ' + n k+q(^ + Wk' - ™k'+q)} 



+ 7]y2 S 4£2 4^ £k+q '' A ~ £ k+q+q',A) (^k+q,A ~ £k+q+q',A) 
qq' 1 1 



X 



{<q',A© q (A,AA)+^ !A e q ,(A,AA)} 



(n) 
M k,A 



and 



2" 

X {™k+q'(™k+q+q' ~ ™k+q ~ D) - m k+q n k+q+q / } 



(ft) 

(n) _ (n) u k,A 



<4a-aa = <q,A + ^© q (A,AA). (65) 

Note that the factor l/2 n was incorporated in f k n qCT , in order to keep the shape of the ansatz 
(|59|) unchanged, 



Cfar W = M k,A C ka + ^ Z^^k^.A ^2 ' a ^)( £ k',A ~ ^k'+q.A) C k , +qQ Ctf/J C k+q7 . 

qk' q a/3-y 

(66) 

For n — > oo, we arrive at the fully renormalized operator 



£(^°°)(A = 0) = W k C k(T + ^ J^kq ^ ^2{#aP ■ £77) fo' - £k'+q) C k , +qQ C k+q7 , 

qk' q a-Pl 



(67) 
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where u k = u k)^o' an d ^k, q = u kq^o- Using H, the expectation values n k and m k as well 
as the spectral function SG(k, uj) can be evaluated. However, due to the strong correlations 
in 7Y, additional approximations will still be necessary. 

To evaluate the spectral function QG(\s.,uj), we start from Eq. ( 1581) for n — ► oo, A = 



3fG?(k,w) 



7?i' 



(68) 



Here c£^* (A — ► 0) is given by Eq. (1671) . The time dependence and the expectation value 
are defined with Ti,, and L is the Liouville operator to 7i. For a state close to half-filling, the 
following relation is approximately valid according to Appendix [B] 



k<7 



Tt, c k(T 



(69) 



It means, in the case that the dynamics is governed by the Hamiltonian TC, in which no 
magnetic interactions are present, a hole can move almost freely through the lattice. Using 
Eqs. ( l67j) and ( 1681) . the spectral function SG(k, uj) then reads 

9fG(k,w) = *£D<?(w-e k ) + 

2 



3L> v ^ 

+ 2iV2 ^ 



qq 



^k,q 

4u>| 



l^k+q' — ^k+q+q' 



X {^k+q+q' m k+q' + ™k+q(-D + ™k+q' ~ ^k+q+q')} 

1 t/q Jq' _, . , _ 

u k,q ^k,q' l £ k+q' — ^k+q+q' ) l^k+q - ^k+q+q' 



(70) 



2 4^4^, 



X {(^k+q' - mk+q)%q+q' ~ ™k+q' (™k+q + D)} 



8{U + £k+q+q' — ^k+q' ~ ^k+q) • 



Note that in deriving Eq. (l70~l) . an additional factorization approximation was used. Thereby, 
an expectation value, formed with six fermion operators, was replaced by a product of three 
two-fermion expectation values. The new quantities nk and rh^ in Eq. (J70|) . 



H 



are again k-dependent occupation numbers for electrons and holes per spin direction, How- 
ever, they are defined with the fully renormalized model H instead of with 7i as in Eqs. (1621) . 
For nk and mk, we use the Gutzwiller approximation 19 



n k = (D-q) + qf(s k ), 
fhk = q(l- /(£k)) with q 



(71) 



1 — n 



8 



1 - n/2 1 - n/2 ' 
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where /(ffk) is the Fermi function, f(ik) = 9(— £k) for T = 0. Note that % is proportional 
to the hole filling 6 = 1 — n. Obviously, the application of c kcr on a Hilbert space vector is 
non-zero only when holes are present. In contrast, does not vanish even at half-filling. 

According to flTUj) . the spectral function 3(j(k, u) consists of two parts: The first one is 
a coherent excitation of energy with the weight u^D. The second part describes three- 
particle excitations. Also note that the sum rule 

/°° n 
dw$G(k,u;) = ([cL,cW]+> — 1 — - — D (72) 

is automatically fulfilled by (!70|) . The sum rule is built in by the construction of the renor- 
malization equations for u^x and v k, q ,A in Appendix O 

For finite temperature, a phenomenological extension of the Gutzwiller approximation 
according to^ will later be used. Here, the Fermi function is replaced by 

f ^ = l+exp^MM)] ' <73) 
where w(k, n) is a weighting function in k-space. It was introduced in22 in order to account 

for an over-completeness in the Gutzwiller approximation. It plays the role of a k-dependent 

effective mass and is a quantity of order 1. 

Finally, note that the static expectation values and m^, defined in Eq. fl62|) . can also 

be evaluated from A(k, uo) or 3G(k,w): 

/OO noo | 

A(k,u)du= - ^ $SG(k,u) du, m k = D - n k . (74) 

-oo J — oo ^ T C 

V. NUMERICAL EVALUATION FOR THE PSEUDOGAP PHASE 

The renormalization equations (jHJ, fl53|) . fl64l) and fl65l) together with (1741 form a closed 
system of equations, which could be solved self-consistently. However, to simplify the numer- 
ical evaluation, we calculate the expectation values in Eq. fj4T|) and Eq. fl53l) with the renor- 
malized Hamiltonian 7i instead of with 7i. Within this approximation and the Gutzwiller 
approximation flTTT) . the renormalization equation for the energy reads 

Jl 



1 . J 

£k,A-AA — £k,A = T7TT7 / TT" ©q(A, AA) (&k+q,A + ^k-q.A ~ ^Vl.x) (S q S_ q ) 

q q> A 



77 ^(2£k',A - ^k'+q,A - £k'-q,A)/k^ L) 



A 

k 



x(£k,A-^k-q,A)VK\ ( 75 ) 
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with 

/a c \ _ 3 ^ 1 VV* s \2 H NL ) f( NL ) 

k' 

Here, ^ s the non-local part of the Fermi distribution, fff^ = 1/(1 + e^ 6 *) — 

(1/N) ^ k 1/(1 + e^ £k ). Remember that the factor g as well as are proportional to 
the hole concentration 8 — 1 — n. Therefore, the renormalization contributions to Eq. f!75|) 
are almost independent of 8 and turn out to be very small. Therefore, from now on, the A 
dependence of and also of lj^x will be neglected. 



A. Zero temperature results 

For the evaluation of the renormalization scheme, we have used a sufficiently large number 
of renormalization cycles in order to obtain self-consistency. We have considered a square 
lattice with N = 40 x 40 sites and a moderate hole doping, such that the system is outside the 
anti-ferromagnetic phase but not yet in the Fermi-liquid phase. Possible superconducting 
solutions are not considered. 

The main feature of the normal state is the appearance of a pseudogap which is ex- 
perimentally observed in ARPES measurements. A small next-nearest neighbor hopping 
if = O.lt and an exchange constant J = 0.2t between nearest neighbors are assumed. The 
inclusion of a non-zero if leads to a Fermi surface (FS), as sketched in the inset of Fig. [3j It 
closely resembles the Fermi surface of non-interacting electrons. The FS is determined from 
the condition = for a fixed value of the electron filling n — 1 — 8. The temperature is 
set equal to T = 0. Let us first concentrate on the u>- dependence of the spectral function 
&Cr(k, a;). In all figures, the symmetrized function will be plotted in order to remove the 
effects of the Fermi function on the spectra. 

Fig. [2] shows the PRM result for SG(k, to) for two different hole concentrations in the 
underdoped regime (a) 8 = 0.03 and (b) 8 = 0.075, for several k-values on the FS between 
the nodal point near (7r/2,7r/2) and the anti-nodal near (tt,0). As the most important 
finding, one recognizes the opening of a pseudogap for both hole concentrations, when one 
proceeds from the nodal towards the anti-nodal direction. On a substantial part of the 
FS, the spectra show a peak-like behavior around uo = 0, indicating a Fermi arc of gapless 
excitations. Note that our analytical results show a remarkable agreement with findings 
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FIG. 2: Symmetrized spectral function 9G(k, u) at T = for two hole fillings (a) 5 = 0.03 and 
(b) 5 = 0.075 along the Fermi surface. The top 9G(k, u) is at the node, whereas the bottom is at 
the anti-node, as defined in the inset of Fig. 

from ARPES experiments in high-temperature superconductors*^^^. Also additional 
peaks are found in the nodal direction at lower binding energies which are enhanced for 
5 = 0.075. In Fig. [3], the pseudogap on the FS is shown as a function of the angle 0, where 
is denned in the inset of Fig. [3j The results are taken from Figs. Et^a) and (b). Note 
that for the smaller hole filling, the length of the Fermi arc becomes smaller, whereas the 
pseudogap becomes larger. This behavior agrees with the known experimental feature of a 
characteristic pseudogap temperature T* which increases with decreasing hole filling^,— . 
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FIG. 3: Pseudogap size A pg from Fig. [2] as function of the Fermi surface angle <f> for 5 = 0.03 
(black) and 5 = 0.075 (red). 

The uj- and k-dependence of SG(k, uj) from Fig. [2] can easily be understood from equation 

(EQD, 

3G(k,u;) = \u k \ 2 D5(u-e k ) 

3D . f / J \ 2 

+ 2iV2 E (4^ l f ' k ' q | 2 ( £k ' ~~ £ k'+q) 2 (^k'+ q m k / + n k+q (£) + fi k / - n k '+ q )) 

qk' L ^ 

H J 5 (w + £k'+ q - £k' - ^k+q) , (76) 

where the dots + • • • indicate additional terms which are less important. First, from the 
renormalization equation fl60l) for w kA , one finds that its original value u 2 . = 1 at A = A is 
reduced by renormalization contributions of order 5~ 2 according to w k A _ AA — w kjA = —a\/S 2 . 
Thus, the weight of the coherent excitation |tt k | 2 becomes small for small 5, so that the 
spectral function SG(k, uj) is dominated by the incoherent excitations in Eq. (I75|) . What 
remains is to show that the different behavior of SG(k, uj) in the nodal and in the anti-nodal 
region can be understood solely from the incoherent part of Eq. (l76j) : 
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First note that the dominant contribution in Eq. (I76j) at small u arises from the small 
q-terms in the sum over q, since in the denominator uj 2 ~ q 2 . In the numerator, the factor 
(e k /— e k / +q ) 2 is also proportional to q 2 , so that the combined prefactor ( J q /4u) 2 ) 2 (£: k /— £ k / +q ) 2 
behaves as ~ q~ 2 . However, the small q terms do not lead to a divergency in Eq. (ITS]) since 
the additional renormalization parameter v 2 . also vanishes for q — > 0. This behavior can be 
verified by a close inspection of the renormalization equations (1601) . (161]) for w^a and Wk,q,A- 
Next, let us use the small q expansion for the energy difference 

£k' - £k'+ q = -2t(q x sin k' x + q y sin k' y ) . (77) 

The excitations from the ^-function in Eq. (1761) are given by 

u = e k / - £ k / +q + e k+q « e k + 2t (sin k x - sin + q y (sin fc y - sin k' y ) } , (78) 

which still depend on k'. There is also a k'-dependent factor in the numerator which con- 
tributes to the intensity, 

Ok' - £k'+ q ) 2 = 4t 2 (q x sin k' x + q y sin k' y ) 2 + 0(q A ). (79) 

Now, we are able to discuss the small tu-behavior of the spectral function 3G(k, w), when 
the wave vector k is varied: 

(i) First, close to the anti-nodal point k = (0, tt), the excitation energy (1751) reduces to 

u = £ k / - e k '+ q + £ k+q « £ k - 2t (g a sin fc^. + g y sin fc£) . (80) 

By comparing Eq. f[8"Uj) with Eq. (1751 . one realizes that the square of the frequency shift in 
Eq. (1801) is identical to the intensity factor (179|) . Thus, excitations with small shifts away 
from the Fermi surface e k = also have small intensities, whereas those with large shifts 
have large intensities. This explains naturally the pseudogap behavior at the anti-nodal 
point, where a lack of intensity is found at u — 0. 

(ii) For the nodal point near k = (tt/2, 7r/2), the excitations have energies 

lo = e k / - e k / +q + £ k+q « £ k + 2i (1 - sin A£) + g y (1 - sin A;^)} , (81) 

whereas the intensity factor is again given by Eq. (1791 . The largest intensity is caused 
by terms in the sum over k' which either belong to the region around k' m (tt/2, tt/2) or 
around k' « (—tt/2, —tt/2). In the first case, the excitations (IHTj) reduces to « e k , whereas 
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the intensity factor (!79l) is given by 4:t 2 (q x + q y ) 2 - Thus, from this k'-region, one obtains 
excitations directly at the Fermi surface. For the second k'-region, the excitation energies 
are given by u ~ £k + 4:t(q x + q y ). The intensity factor is the same as before. Thus, similar to 
the anti- nodal point, the square of the excitation shift away from the Fermi surface = 
is proportional to the corresponding intensity. Therefore, from these k'-terms no intensity 
is expected at u = 0. To summarize, an excitation peak at uo = is expected for wave 
vectors k at the anti- nodal point from the first k'-regime, discussed above. In contrast, 
for wave vector k at the anti-nodal point a pseudogap arises. This explains the pseudogap 
behavior of the ARPES spectral function and leads to an understanding of the spectra of 
Fig. [2J In Fig. HJ the spectral function is plotted for a larger hole concentration 5 = 0.09. 
The remarkable new feature is the occurrence of a narrow coherent excitation at ou — 0. 
Note that for this hole concentration, the weight _D|-Sk| 2 of the coherent excitation is no 
longer negligible as in the preceding cases since the renormalization contributions ~ 1/S 2 to 
u 2 , A are less important for larger 5. By increasing 5, the coherent peak gains weight at the 
expense of the incoherent excitations. We also expect a broadening of the coherent peak 
due to a coupling to other degrees of freedom such as phonons or impurities. 

In Figs. |5]^a) and (b), the spectral functions are shown for two different cuts in the 
Brillouin zone. In both figures, k x is fixed and k y is varied thereby crossing the FS. In 
panel (a), where k x = tc, the cut runs along the anti-nodal region through the FS at kp « 
(7r,0.077r). Note that the pseudogap is restricted to a small k-range around the anti-nodal 
point. It disappears for larger k y values away from the anti- nodal point, in agreement with 
the earlier discussion on the origin of the pseudogap. The spectra along a cut in the nodal 
region are shown in panel (b), where k x = tt/2. Apart from the dominant excitation which 
corresponds to the gapless excitation on the FS in Fig. [2j also weaker excitations are found at 
lower binding energies. The complete peak structure is shifted almost unchanged through 
the FS, when k y is varied. The energy distance between the primary and the secondary 
peak slowly decreases by proceeding along the FS from the nodal point to the anti-nodal 
point, until finally both peaks disappear when the anti- nodal region is reached. Such a 
double-peak structure with the same properties along the FS was observed in underdoped 
cuprate superconductors^. Finally, one point might still be worth mentioning. For fixed 
u, the spectrum in k-space is much broader than what one would expect for free electrons. 
Thus, the electron occupation (cj^Ckcr) = / du(l + e^)~ 1 Q ; G(k, lu) depends only weakly on 
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FIG. 4: Same quantity as in Figs. Gla) or (b) for a larger hole doping of 5 = 0.09. 

k. This feature is consistent with the former expression (1711) for n^, where the Gutzwiller 
approximation was used. Remember that the expectation value was defined with the 
renormalized Hamiltonian 7i and not with Ti. 

B. Finite temperature results 

Next, we discuss the influence of the temperature on the one-particle spectra in the 
normal state. For the hopping to next nearest neighbors, we use a somewhat larger value 
t' = 0.41 This leads to an enhanced curvature of the Fermi surface, as it is observed 
in most of the copper oxides superconductors. The other parameters remain unchanged. 
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FIG. 5: Spectral functions 3 : G(k, oS) for two fixed k x values: (a) k x = n and (b) k x = ir/2 and 
different values of k y , thereby crossing the Fermi surface. The hole filling 5 = 0.075 is the same as 
in Fig.E^b). 

Fig. [6] shows the symmetrized spectral function SG(k, u) for three different temperatures 
(a) T = 0, (b) T = 0.04t, and (c) T = 0.08i. The hole concentration for all curves is 5 = 0.04. 
Possible superconducting solutions are again suppressed. The results are shown for different 
k- vectors on the Fermi surface between the nodal (top) and the anti- nodal point (below). 
For all temperatures, a separation of the Fermi surface into two segments is found, as it 
was already discussed in the foregoing section: (i) For k-vectors around the nodal points, 
3G(k, uj) shows strong excitations at uj = (black curves). They form the Fermi arc. (ii) The 
other segment is given by k-vectors, for which 3G(k, uj) shows a pseudogap around uj = 
(red curves). From Figs. E]^a)-(c), one can see that the length of the Fermi arc increases with 
increasing temperature. This increase is equivalent to a reduction of the pseudogap region. 
For instance, for the largest temperature T = 0.08t, the pseudogap is restriced to a quite 



29 
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FIG. 6: Symmetrized spectral function 3?G(k, uj) at doping 5 = 0.04. for three temperatures (a) 
T = 0, (b) T = 0.04t and (c) T = 0.08i for k- values along the FS. The other parameters are 
t' = OAt and J = 0.2t. The top 9G(k, uS) is at the node, whereas the bottom is at the anti-node. 
A possible superconducting solution was suppressed. 

small region around the anti- nodal point. Note that this temperature behavior is in good 
agreement with recent ARPES experiments 9 . A comparison of the spectral functions at 
the anti- nodal point for three different temperatures (lowest curves in Figs. Efa)-(c)) shows 
the influence of T on the pseudogap: With increasing T, the pseudogap is filled up with 
additional spectral weight, whereas the magnitude of the gap (i.e. the distance between the 
maxima on the cu-axis) remains almost constant. Also this temperature behavior is verified 
experimentally 9 . A characteristic temperature T* can be defined at which the pseudogap 
is completely filled up, and the Fermi arc extends over the whole Fermi surface. This 
temperature T* was already introduced above and is called pseudogap temperature. For the 
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FIG. 7: Pseudogap A pg from Fig. [6] as a function of the Fermi surface angle (f> for T = (blue), 
T = 0.04t (red), and T = 0.08t (black). 

present case, T* is approximately T* pa O.lt. 

The pseudogaps, taken over from Figs. [B^a)-(c), are shown in Fig. [7| for three different 
temperatures as function of the Fermi surface angle 0. Note the strong increase of the 
pseudogap at a finite Fermi angle which depends on the temperature. This particular angle 
marks the transition between the Fermi arc and the pseudogap section. At T = 0, it is about 
25 degrees and moves towards the anti-nodal point for higher temperatures. From Fig. [3, 
one may also deduce that the length of the Fermi arc approximately increases linearly with 
T. Also this feature is consistent with ARPES experiments^. 

To discuss the influence of S on the temperature dependence, in Fig. [8] the symmetrized 
spectral function SG(k, u) is shown as function of u for two different temperatures T = 
(black) and T = 0.08t (red) and for five different hole concentrations between 5 = 0.04 
(bottom) and 5 = 0.075 (top). The k- vector is fixed to the anti-nodal point on the FS. The 
curves for T = (black) show a decrease of the pseudogap with increasing hole concentration 
until it vanishes at S ~ 0.075. For the higher temperature T = 0.08t (red), the pseudogap 
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FIG. 8: Symmetrized spectral function 9G(k, u) for fixed k value on the anti-nodal point for five 
different hole concentrations from 5 = 0.04 (bottom curves) to 5 = 0.075 (top curves). In each 
case, the results are shown for two different temperatures T = (black) and T = 0.08i (red). For 
the coherent excitations ~ |"Uk| 2 , the same broadening has been taken for each 5- value. 

vanishes already at a lower hole concentration of 5 ~ 0.06. This verifies the experimentally 
known decrease of the pseudogap temperature T* with increasing hole concentration. 

The doping and temperature behavior of 3G(k, u) can be understood on the basis of the 
former result (1761 for the spectral function. First, in Fig. |9l the parameter is shown as a 
function of 5 which shows a strong increase with the hole concentration. According to the 
first line in Eq. (1761 . agrees with the amplitude of the coherent excitation. Therefore, in 
Fig. [6] for instance, the weight of the coherent excitation ~ |Mk| 2 is negligibly small for the 
smallest hole concentration 5 = 0.04, and the spectrum is dominated by the incoherent part 



32 



0.8 - 
0.6- 

IS 

0.4- 
0.2 - 




n i i i i i i i i i i i 

0.05 0.1 0.15 0.2 0.25 

hole doping 

FIG. 9: The renormalized amplitude of the coherent excitation in Eq. (|76p is shown as a 
function of the hole concentration 5. The k- vector is fixed to (0, 7r). 

of Eq. (1761) . In contrast, for sufficiently large 5, a coherent excitation at u = is expected, 
when k is fixed to the Fermi surface. This behavior is for instance realized in Fig. HI Note 
that an additional broadening of the coherent excitation should be included, which follows 
from the scattering of the charge carriers at additional phonons or impurities. In Fig. [HI 
this broadening was assumed to be T-independent and was set equal to O.lt. Therefore, 
the following doping behavior can be deduced from Fig. [HJ For small hole concentrations 
5 (5 <C 0.07), the spectrum at T = is dominated by the incoherent excitations with a 
pronounced pseudogap around the anti-nodal point. For intermediate hole doping (5 ~ 0.07), 
the spectrum is a superposition of a coherent and of incoherent excitations. Both parts are 
of the same order of magnitude for an intermediate doping. The incoherent part has still a 
pseudogap which is partly compensated by the broadening of the coherent excitation. For 
larger doping 5 > 0.07, the spectrum mainly consists of a coherent excitation around uj = 0. 
With respect to temperature, the coherent excitation is almost unaffected by T, whereas the 
pseudogap is filled up due to the temperature-dependent shift of the Fermi surface, as will 
be explained below. 
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To understand the T-behavior of the spectral function, keep in mind that and therefore 
the weight of the coherent excitation in QfC? (k, a>) , is almost independent of T. Moreover, 
the total spectral weight, to which coherent and incoherent excitations contribute, is T 
independent. This follows from the sum rule (1721) . since the total electron number is fixed. 
Thus, except of minor changes, the overall temperature dependence of 3G(k, uj) is expected 
to be weak. Instead, the main reason for the T-dependence can be traced back to a change 
of the Fermi surface with temperature. Consider a k-vector on the Fermi surface at the 
anti-nodal point, kp = (tt, k y F ), where the x-component is fixed to kp = ir. By varying the 
temperature, one finds that the magnitude of the y-component kp increases almost linearly 
with T. Due to this shift of the Fermi energy with T, also the positions of the incoherent 
excitations at to = are shifted. In this way, one understands that the pseudogap is less 
pronounced for higher temperatures, when hp is fixed to kp = n. A similar behavior of the 
pseudogap was found before in Fig. [5j There, the spectral function is shown for fixed k x = n 
and different values of k y , when the temperature is fixed. Also in this case, the pseudogap 
is suppressed for larger values of k y . Finally, note that kp also strongly depends on the 
nearest-neighbor hopping t' . For small t' , the pseudogap is more pronounced than for larger 
values of t'. This can be seen by comparing the spectrum in Fig. [2] (with t' = O.lt) with 
that of Fig. El where t' = OAt. 



VI. CONCLUSIONS 



In this paper, we have given a microscopic approach to the pseudogap phase in cuprate 
systems at moderate hole doping. Thereby, a recently developed projector-based renormal- 
ization method (PRM) was applied to the t-J model. The pseudogap, which is found in 
ARPES experiments, can be traced back to incoherent excitations in the one-particle Green 
function. It can neither be explained by a competing order nor as a precursor of super- 
conductivity. Instead, the pseudogap phase is an intrinsic property of the cuprates close to 
half-filling. In a subsequent paper—, we shall show that a transition to a superconducting 
phase occurs in the formalism either by lowering the temperature or by approaching an 
appropriate doping range. 
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APPENDIX A: DERIVATION OF THE SPIN SUSCEPTIBILITY x (q, w) 

The derivation of the spin susceptibility x(q, uj) in Eq. (1241) for the system, described 
by the Hamiltonian Ti^ — T~L t + 7ij , is based on the Mori-Zwanzig projection formalism. 
This formalism allows to derive exact equations of motion for an appropriately chosen set 
of relevant operator variables {A*}, 

jA a {t) = iJ2Mt)^a- I J2Mt-t')Zpa(t'W + F a (t), (Al) 
P ^° P 

where the dynamics of the set A a (t) should be governed by Ho, i.e. A a {t) = e^ Wo< A a e~^ not . 
The quantities iu) a p, S Q/3 (t), and F a (t) are called frequency matrix, selfenergy, and random 
force 

in aP = J2x^(A a \Ap), = 5>S (i 7 |Qe iQLoQ *QM (A2) 

7 7 

F a {t) = i e iQL ° Q *QL A a . 

Here, A a is the time derivative of A a , defined by A a = i L A a , and Xap ls the inverse of the 
susceptibility matrix x a p = (A* I As)- In Eqs. ( \A2\i , we have also introduced a scalar product 
between operator quantities A and B, 

rP 

(A\B) = / d\(A j e- XLo B) , (A3) 
Jo 

where the expectation value (• • ■ )o is formed with 7i and L is the Liouville operator, which 
corresponds to 7io- In ^ap(t) the quantity Q is a projection operator which projects on the 
subspace of all operator variables which are 'perpendicular' to the set {A*}, i.e. 

Q = l-^K)x^l- (A4) 

a 

To use the general projection formalism to derive x(q> u), we have to choose an appropri- 
ate set of relevant operator {AA In our case, this set is given by S q and its time derivative 
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S q , i.e. 



{AJ = {S q ,S q }. (A5) 
From the equations flAlj) . one easily derives the following two equations: 

jS^t) = S q (t), (A6) 
^S q (i) = - W ;S q (t)-^t'S q (f-t')E q (0+F q (f), 
where the frequency and the selfenergy in the second equation are given by 

< = W^Tv E q (t) = ^^(SjQe^^QSq) (A7) 



and the random force is F q (t) = e jQL ° Q *QS q . The projector Q projects perpendicular to S q 
and S q . In deriving the equations (1A6I) . we have also used (S^ \S£) = i([S" , S£\)q = (for 
all u, fj, = x,y,z), which follows from the exact relation (A|L S) = ([z^, .B]) . To find the 
dynamical susceptibility x(q, u), we multiply both equations flA6j) with the 'bra' (S q | and 
go over to the Laplace transform. Using (S q |F q ) = 0, we obtain 

— oj 2 

= 2 q v / x X(q)- (A8) 

Here, x(q) = O^ql^q) i s the static spin susceptibility and £ q (w) is the Laplace transformed 



selfenergy 



E < H = (S^ (§JQ ^QUQ-„, Q ^' (A9 » 



To proceed, we have to evaluate the second time derivative S q 
S q = ~ £ t\ (e^ - e^) (5> (.) - &> o (0) 

52 ^^(e^-e^) (A10) 

x (ct a P a (Z) c i/3 + c^-a 5 ?^) + [c\ p V a {l) C ja + clpS^Cm^ } , 



where only the dominant part of the hopping Hamiltonian TCt was taken into account. The 
first term on the right hand side of Eq. (1A10I) enters from a twofold hopping to a neighboring 
site and back. By replacing the two projectors Vo(i) and Vq(1) by their expectation values, 
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we come back to the former equation of motion (122]) . Therefore, we can conclude that the 
frequency term u; 2 , defined in Eq. ( 1 A7I) . agrees with the former frequency term o> 2 from 
Eq. 



= = 2P (*J =0 - tj) > 0. (All) 

The second contribution in Eq. (lAlOf) describes a twofold hopping away from the starting 
site and agrees with the quantity QS q in the selfenergy, 

qs q = ( eiqRi - e " Ri ) x ( A12 ) 

x Q {a Q/3 (cj a £> a (7) cm + c)_ a S?c if ^ + a*^ (c\ V a (l) c ja + c\ Sf a c TO _ a J 



In order to obtain a rough estimate for the selfenergy £ q (a>), we neglect the spin flip operators 
in Eq. (1A12I) and replace the local projectors T> a (i) and T> a (l) as before by their expectation 
value D. By introducing Fourier transformed quantities, we find 

Q ^ q = 27W ^ 9aP ((£k+q " £fc)2 " 2(tq=0 " tq) ^ Q (A13) 

a/3 



The selfenergy then reads 

D 2 



^EEE^-^' ( A14 ) 



q y ^ c i q \ 4iV ' — " ' — * - — •• 

(^D q |D q J kk' a/3 a'/3' 

x [(e k+q -e k ) 2 - 2(t 2 =0 - Q] [(e k , +q - eu) 2 - 2(t q=0 - t\)} 



x (4 +q , a ^i ^ _ qi^q _ iv Q 4 +q ,^, 



In the final step, we factorize the two-particle correlation function in Eq. (IA14I) in a product 
of one-particle Green functions. A straightforward calculation leads for the imaginary part 
of the selfenergy to 

9E q ( w ) = _^!_A^[ (£k+q _ £k) 2_2(t q=0 -t q )] 2 SAf k (q, W ), (A15) 

(SqlSq) k 



SM k (q,cu) = / duo 



(3oj vr 7_ 00 1 + e -/J(«+a) 1 + 

Here, SG^^u;) is the imaginary part of the one-particle Green function, formed with the 
Hamiltonian Tio, 

POO 

Gf M = i / dt([c Ka (t),clj + ) c -*(«-*»)*. (A16) 
Jo 
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Finally, we have to evaluate the denominator (SqISq) of E q (o;). Proceeding in analogy to 
the evaluation of £ q (o>), we find 

3 . 

(Sq|Sq) = ^Xl( £k + t l~ £k ) 2 ( a k+q, a Ck/3|cj c + q , a Ck ( 0) (A17) 
k 

with 



APPENDIX B: FACTORIZATION APPROXIMATION FOR S n A S 



q,A>->-q,A 

The aim of this appendix is to simplify the operator product Sq^S^qA in the expressions 
for Ho,a and TC\ t \ from Sec. IIV A[ 

^0,A = ^ "77 ( S q ■ S_ q H S q A ' S_ q A J , 

q 2 V ^q,A J 

7~tl,\ = ( Sq ■ S_ q Sq ; A • S_q 5 A J • 

q 2 V ^q,A / 

This will be done by use of a factorization approximation. Using for the time derivative 

we first can rewrite S qi AS_ qi A as 

Sq,AS_ q ,A = ^ E • £ / ^ A«>' <,H - c **0 x 

x E w(^ iqR! - e- iqRm ) 4^4^- ( B1 ) 

Using a factorization approximation, the four-fermion operator on the right hand side can 
be reduced to operators cj^c^ which will lead to a renormalization of e^. Thereby, we have 
to pay attention to the fact that the averaged spin operator vanishes ((Sj) = 0) outside the 
antiferromagnetic regime. Moreover, all local indices in the four-fermion term of Eq. ( 1B1I) 
should be different from each other. This follows from the former decomposition of the 
exchange interaction into eigenmodes of L t in Sec. IIV A\ where we have implicitly assumed 
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that the operators S q .A and S_ q .A do not overlap in the local space. Otherwise, the decom- 
position would be much more involved. However, it can be shown that these 'interference' 
terms only make a minor impact on the results. For the factorization, we find 



X 



^2((^€i/3)nl) (c\ a ci a ) NL + ^2((c\ a ci a ) NL ) (c jl3 c m p) NL \ , (B2) 

2 /3 J 

where we have neglected an additional c-number quantity, which enters in the factorization. 
The attached subscript in (■ ■ - )nl on the right hand side indicates that the local sites of the 
operators inside the brackets are different from each other. Note that sums over spin indices 
in Eq. ( 1B1I) have already been carried out. Fourier transforming Eq. ( 1B2I) leads to 

Sq,AS_ qi A = -— ^( £k ' A ~ £ k-q,A) 2 ((4-qaCk-qa)ArL) (c ko .C k <x) NL, (B3) 
kcr 

where we have defined 



N 
k' 



Using Eq. flB3j) together with Eq. (l5Tj) . one is led to the renormalization result fl53|) of e 
to first order in J. 

In the following, let us simplify the notation and suppress the index A in S q A, £k,\, and 
also in cj q .A- With this convention, we shall use the factorization (1B3P in order to derive the 
renormalization (HIT) for e^x in second order in J. We start from expression (flOl) for the 
renormalized Hamiltonian T^a-aa ^ n secon d order 

WA-AA = {©(A-|2^ q ,A|)-i} [XA,AA,AA(q)+4,A(q)] + E J ^ X ^A,A,A] 



£ J q q (A, AA) (j[X XAX , S q • S_ q ] + ^[Xa,aa , S q • S-Sj 



(B4) 



where in the first line we have already used P^a.aa , T~tt,\) — — J2q @q(^; AA) (•^•i.aCq) + 
*4i a(^))- Next, we have to evaluate the commutators of Xa,aa with S q • S_ q and S q • S_ q . 
Using [S^ q , Sq] = ^ E q(T ( 2£: k - £k+ q - e k _ q ) c k <A<7, (f = 2/, 2;), and Eq. ([39]), we find 



[X A ,AA , S q • S_ q ] = A 0q(A)AA) ( 1 ^(2 £k - £k+ q- £k _q)(4 CT C kCT ) 

q \ ka 

+ x4 e q( A > AA ) ( s i ' s -q) 4 5Z( 2£k ~ £k +i ~ £k -i) a ^ akCT 



S q ■ S_ q 



1 kcr 
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[X XAX , S q • S_ q ] = -A 0q(A;AA ) ^l^( 2£k -e k+q - £k _ q )(cj C(T c k(T )j S q -S_ 



J, 



5L 6 q (A, AA) (S q • S_ q ) 1 ^(2e k - e k+q - £k _ q ) c k(T c k(T . (B5) 



1 " ka 

Note that in (IB5I) already a factorization approximation was used. With the relations (1B4|) 
and (IB5I). we obtain 



^ J^(2e k - £ k+q - £ k - q )(c kcr Cka) 

kcr 



Sq ' S q 



+ (S q • S_ q ) — ^2(2e k - £ k+q - £ k _ q ) c[ a C ka J 
kcr / 



(B6) 



J. 



^(2e k - e k+q - £ k -q)(c kCT c kCT ) 

kcr 



Sq ' S_ q 



+ (S q ' S_ q ) — ^(2£ k - £ k+q - £ k 
kcr 



-qj ^kcr^kcr 



In a final step, we factorize ~ S q ■ S_ q according to (|B3 



n 



(2) 

A— AA 



3£ 



^e q (A,AA) 



- ^(2£k - ^k+q - £k-q)(c k crCkcr) 

kcr 

+ ( S q ■ S -q) Jr 5Z( 2£ k - ^k+q - £k-q) 4cAcr J 

kcr / 

E( ^ A ' AA ) ^ ' S-q) ^ Z>k - ^+q " £ k -q) C^Ck 



Sq " S_ q 



q q 



kcr 



+ 4B4) 2 e,(A,AA) 



qcr 



X^ 2£k ' _ £k ' + C l _ £ k'-q)(C k 'cr'Ck'cr') 



X 



X E( £k ~~ £ k-q) 2 ((c k -qaCk-qa)7VL) (4cA<xW 



(B7) 



kcr 



From (IB7I) . the renormalization equatiuon (HTj) for £k,A-AA can immediately be deduced. 

APPENDIX C: RENORMALIZATION EQUATIONS FOR FERMION OPERA- 
TORS 

The aim of this appendix is to derive the renormalization equation for the fermion op- 
erator Ckcr(A) = e Xx dk a e~ Xx in second order in J q . As before, we shall suppress the index 
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A everywhere in S qj A, ^q,A, and in order to simplify the notation. Let us start from 
an ansatz for c kcr (A) after all excitations with transition energies larger than A have been 
integrated out. It reads 

c k<T (A) = u k>A c kCT - i y2 ©(|2o)q| - A) v k ^ x [S q • S_ q + S_ q • S q , c ka ]. (CI) 

q 1 

In Eq. (ICip . the parameters w kj A and f k , q ,A account for the A-dependence. Note that the 
operator structure in Eq. (IClj) corresponds to that of the first order expansion for c k(T (A) ~ 
c k(J + [^A,c ko -]. Here, X\ has the same operator form as the generator X Aj aa in Eq- (1381) . 
Due to construction, the q-sum in Eq. (101 j) only runs over q- values with excitation energies 
|2cu q | larger than A. This is assured by the B-function in Eq. (ICip . For simplicity, in the 
following we agree upon to incorporate the 0-function in u k , q ,A- Thus, we can write 

C ko -(A) = M kiA C k(J - Zy^ k ,q,AVT^ (([S q ,C k(J ] • S_ q + S_ q • [S q , C k(T ]) 

q q 

+ (S q - [S_ q ,c k(J ] + [S_ q ,c kCT ] ■ S q )). (C2) 
For the additional renormalization from A to the reduced cutoff A — AA, we have 

c k(J (A - AA) = e x ^ c k(T ( A) e~ x ^ = (C3) 

= Uk^^c^e-**^ - i V ^ k , q ,A A e x ^[S q ■ S_ q + S q ■ S q , c ka }e~ x ^\ 

q AuJ * 

where Xa,aa is the generator from Eq. ( |38l) . 

X XA x = -^i|- e ^ A ' AA ) (SqS_ q + S q S_ q ). 
q q 

First, let us expand the term ~ u^x in Eq. (1C3I) . 

e XA ' AA c kCT e~ x ^ = c k(T + [X a ,aa , c kCT ] + -[Xa,aa, [X a ,aa, c k(7 ]] + ■ ■ ■ . (C4) 
Here, we can combine the second term in Eq. (1C3I) with the second part in Eq. (1C2I) . 
c k(J (A-AA) = (« klA + ---)c k «x (C5) 

A 



- i y^(^ k ,q,A + M k ,A @q(A, AA) H ) ^ [S q ■ S_ q + S_ q ■ S q , C k(T ] + 

q 1 



where the dots (+•••) mean additional contributions from higher order commutators with 
^a,aa- On the other hand, c k(T (A — AA) should have the same form as the ansatz (1C1I) . when 
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A is replaced by A — A A, 

c k(T (A - AA) = u k ^-AxCka - ^y^k,g,A-AAj^ [S q ■ S_ q + S_ q ■ S q , c ko ]. (C6) 

q q 

The comparison of Eqs. (1C6|) and (IC5I) immediately leads to the renormalization equation 
([EED for Uk, q ,A, 

^k, q ,A-AA = ^k, q ,A + W k) A @ q (A, AA), (C7) 

where we have restricted ourselves to the lowest order contributions in X\,ax- Furthermore, 
we have exploited the very weak A-dependency of £k,A and & q ,A- 

The renormalization equation for the second parameter requires the evaluation of 
higher order commutators in Eq. (1031) . Alternatively, we can start from the anti-commutator 
relation (j3J) 

[C(A)A.(A)] + = ^Xy*2>,(i)e-*» = ^X>"W 

i i 

with Po-(i) = 1 — Tii-o-, where in the last relation [Xx,'^2 i 'Dff(i)] = was used. When we 
take the average, we obtain 

= (^(i)) =: (C8) 

In order to evaluate the anti-commutator in Eq. (1C8|) . we have to insert the former ansatz 
(IC2p for Cko-(A). Here, we make an additional approximation by taking into account only 
the two first terms in Eq. ( 1C2I) . The remaining terms have explicit spin operators S q . In 
the commutator of Eq. (1C8jl . they lead to additional contributions with one or two spin 
operators. Outside the antiferromagnetic phase, no magnetic order is present and also spin 
correlations are weak. Therefore, it seems reasonable to neglect these terms. Thus, we can 
approximate c ko -(A) by 

Ck CT (A) = u KX c ka - i ^2 ^k,q,A^- ([S q , c k(T ] ■ S_ q + S_ q ■ [S q , c k(T ] ) (C9) 

q q 

= U KX Cka + X^ k >1' A 4£2 X^"' 9 ' X^' ~ £k '+q) 4'+qa £ k'/3 C k+q7 . 

q 1 a/37 k' 
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Inserting Eq. (IC9I) and c ko .(A) into Eq. (IC8I) . we obtain 

q'q q ' q a' ,0' ,j' a,P,-y 

X (£k" - ^k"+q') (£k' - ^k'+q) ( [Ck+q'Y^k"^' ^k"+q'a' , C^'+qa ^k',3 C k+q7 ] + ). 

k',k" 

(CIO) 

To find the renormalization equation for u^x-ax, we use the same equation, thereby replacing 
A by A — AA. We then obtain 

D = I^A-AAi^+T^EKq'-A +<A@q'(A,AA)) K,q,A+U k , A q (A,AA))-^-A 

v > q'q q' q 

X ^ ^(^W ■ Vj'*){Vai3 ' <?a 7 ) (Cll) 

a',/3',7' a,/3,7 

X 

k'.k" 



^ ( £ k" - ^k"+q') Ok' - £k'+q) ( [c k +q' 7 'C k "/3' Ck"+q'a' , 4'+qa £ k'/3 C k+q7 ] + ) 



where we have inserted the former renormalization result ( 1C7I) for fk, q ,A-AA- Restricting 
ourselves to the lowest order contributions in J q , we can subtract Eq. (1C10[) from Eq. (101 ip 
and obtain the renormalization equation which connects w^a-aa with w kj A, 

|Mk,A-AA| 2 ^ = |Wk,A| 2 £ - Yl Yl Yl (^'«' ' ^Ofotf ■ ^7) ( C12 ) 

q'q q ' q a',/3',7' «,/3,7 

X {|M k ,A| 2 ©q'(A, AA)6 q (A, AA) 

+ («k,A< q ' j A e q( A > AA )+<A^k,q,A©q'(A,AA))} 

X ^2 ( £k " - ^k"+q') (^k' - ^k'+q) ( [4+q' 7 4"/3' C k '' + q' a ' , C k , +qQ C k+q7 ] + ). 

k',k" 



What remains is to evaluate the commutator in Eq. flC12j) . In a final factorization approxi- 
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mation, we find 

I 1 2 I 1 2 1 \ / \ 2 \ I -> -> 1 2 

q q a,P,l 
X0 q (A,AA) {|Mk,A| 2 + K,A< q , A +Mk,A^k,q,A)} 

X E^' ~ £k '+q) 2 { ™k+q(™k' + -D) + ^k'+q("ik' - ™k+q) } 
k' 

+ (2^52 E ^4 E • • ^) 

V / q/ q q' q q,^^ 

x{|n kiA | 2 0q,(A,AA)0q(A,AA) 

+ Ka <,q', A e q (A, AA) + u* KX ^ k ,q, A q ,(A, AA)) } 

x (^k+q — ^k+q+q') (^k+q' _ ^k+q'+q) {^k+q' (fik+q + D) 

+n k+q+q /(m k+q - n k+q /)}. (C13) 

Summing over the spin indices and exploiting that and fk,q,A are rea l) we arrive at 
expression (j50|l . 
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